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Abstract 

NPL  has  developed  a  Kalman-filter-based  clock  algorithm  for  combining 
measurements  from  its  three  active  hydrogen  masers.  The  algorithm  is 
designed  to  produce  a  near  optimal  composite  when  the  dominant  noise 
process  present  is  flicker  frequency  modulation  (FFM).  The  FFM  is  modelled 
approximately  by  a  linear  combination  of  Markov  noise  processes.  Each  Markov  process  is 
included  in  the  Kalman  filter  and  contributes  an  additional  component  to  its  state  vector. 
Both  the  validity  of  the  model  and  the  effectiveness  of  adding  these  additional  components 
to  the  state  vector  are  examined.  The  performance  of  the  new  algorithm  is  examined  when 
applied  to  simulated  measurements  and  also  to  measurements  obtained  from  NPL’s 
hydrogen  masers. 


1.  INTRODUCTION 

The  National  Physical  Laboratory  (NPL)  maintains  and  develops  the  UK’s  national  time  scale 
UTC  (NPL).  NPL  has  a  requirement  to  both  improve  the  stability  of  UTC  (NPL)  and  to  steer  UTC 
(NPL)  as  close  to  UTC  as  possible.  The  improved  stability  of  UTC  (NPL)  is  required  for  the 
evaluation  of  new  caesium  fountain  clocks  [1]  and  optical  frequency  standards  at  present  under 
development. 

Clock-combining  algorithms  are  of  two  types:  those  based  on  the  use  of  a  Kalman  fdter,  for  example 
the  GPS  composite  clock  algorithm  [2],  and  those  based  on  a  weighted  average  of  the  clock  data,  for 
example  the  NIST  ATI  algorithm  [3].  Most  operational  clock  algorithms  used  in  primary  timing 
laboratories  are  of  the  latter  type.  There  have  been  problems  in  the  implementation  of  Kalman-filter- 
based  clock  algorithms.  This  is  because  the  physical  parameters  estimated  by  the  state  vector  are  only 
partly  observable  from  the  measurements  made  between  individual  clock  pairs.  Greenhall  [4]  has 
recently  overcome  these  problems  through  the  development  of  a  method  described  as  covariance 
reduction.  The  resulting  clock  algorithm  has  the  potential  to  produce  a  composite  clock  with  a 
stability  that  is  close  to  optimum  for  all  averaging  times. 

Kalman-fdter-based  clock  algorithms  model  the  clock  noise  as  a  linear  combination  of  White 
Frequency  Modulation  (WFM),  Random  Walk  Frequency  Modulation  (RWFM),  and  Random  Run 
Frequency  Modulation  (RRFM).  By  applying  the  “three-cornered-hat”  technique  to  NPL’s  three 
active  hydrogen  masers,  it  is  possible  to  obtain  plots  of  Logio(cTy)  against  Logio(x)  for  each  of  NPL’s 
masers.  The  results  are  shown  in  Figure  1.  The  plots  reveal  a  complex  noise  structure,  where  some  of 
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the  noise  present  may  be  better  described  as  Flicker  Frequency  Modulation  (FFM),  for  which  the 
curves  have  an  underlying  gradient  of  zero,  in  addition  to  WFM  and  RWFM,  where  the  underlying 
gradients  are  — 1/2  and  1/2 ,  respectively. 
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Figure  1.  Plots  of  Logio(Cfy)  against  Logio(x)  obtained  from  measurements  made 
between  pairs  of  NPL’s  masers  and  for  the  individual  clocks  estimated  using  a  three- 
cornered  hat. 


In  this  paper,  we  extend  the  work  of  Greenhall  [4]  to  develop  a  clock  algorithm  that  aims  to  provide 
close  to  optimal  perfonnance  in  the  presence  of  FFM.  The  FFM  is  modelled  as  a  linear  combination 
of  Markov  noise  processes.  Each  Markov  process  is  included  in  the  Kalman  filter  and  contributes  an 
additional  component  to  the  state  vector.  The  properties  of  the  resulting  clock  algorithm  are 
examined  using  both  simulated  and  real  clock  data.  An  iterative  method  has  previously  been 
developed  to  determine  the  magnitude  of  noise  parameters  from  a  measurement  time  series  [5].  This 
method  is  extended  and  applied  to  a  “three-cornered  hat”  of  hydrogen  maser  measurements  to 
estimate  the  required  noise  parameters. 


2.  MODELLING  FLICKER  FREQUENCY  MODULATION  AS  A 
LINEAR  COMBINATION  OF  MARKOV  NOISE  PROCESSES 

Flicker  Frequency  Modulation  (FFM)  may  be  considered  to  be  a  power-law  noise  process  for  which 

2  I  |3  2 

the  spectral  density  function  of  the  phase  measurements  varies  as  S^{f)  =  O'  l\^7lf\  ,  where  o  is 

the  variance  of  an  underlying  white  noise  process  and  f  is  the  Fourier  frequency  [6].  Recently,  FFM 
has  been  described  in  term  of  fractionally  integrated  noise  processes  [6,  7] .  In  this  work,  we  adopt  an 
older  approach  initially  studied  by  Mandelbrot  [8],  where  FFM  is  described  as  a  linear  combination  of 
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Markov  noise  processes.  Let  {x,.  /  =  1,  ...,m}  be  a  time  series  of  measurements  ofx  made  at  uniformly 
spaced  times  i  =  i,...,m}.  A  Markov  noise  process  takes  the  following  form: 

XM=kx,+e,,  (1) 


where  k  is  a  constant  for  an  individual  Markov  process,  6^,6 j  are  independent  random  errors,  such 
that  cov  (gj. ,  gy  )  =  0,  i  ^  j  ,  with  normal  distribution,  zero  expectation  E  (e)  =  0,  and  variance  (T^ . 
The  variance  of  x,  as  i  approaches  0°  is  given  by  U,  where: 


U  = 


{\-ey 


(2) 


Assuming  k'^  <1  then  the  autocorrelation  function  p{l)  is  given  by: 


p(l)  =  exp(-f?  /) 


(3) 


where  /  is  the  lag  and  k  =  exp(— f?ro)  ,  To  is  the  time  interval  between  successive  points.  In  the  work 
undertaken  here,  FFM  is  modelled  as  a  linear  combination  of  h  Markov  processes,  where: 

h 

yfa=Hyj,.  ,  and  f/,- =  var(j^.)  (4) 

M 

yjy  is  the  i*  element  of  the  7*  Markov  process,  and  yQ  is  the  i*  element  of  the  resulting  approximation 
to  FFM  frequency  measurements.  We  choose  the  magnitude  of  the  Markov  components  so  that  the 
variances  of  the  component  Uj  =  var  (yj)  are  equal.  This  will  result  in  time  series  where  plots  of 
Logio(cty)  against  Logio(x)  will  have  a  gradient  of  approximately  zero. 

The  properties  of  FFM  modelled  in  terms  of  a  linear  combination  of  Markov  processes  are  shown  in 
Figures  2,  3,  and  4  below.  Four  Markov  processes  are  used  to  model  FFM  with  =  Rj  /S  and 

R^  —  0.75,  Tg  =  1.  Curves  of  a  typical  FFM  simulation  and  their  component  Markov  processes  are 

shown  in  Figure  2.  Plots  of  Logio(ay)  against  Logio(x)  are  shown  in  Figure  3,  based  on  theory 
(continuous  lines)  and  simulated  data  (marks),  and  including  results  for  the  component  Markov 
processes. 


When  simulated  data  are  used,  the  initial  values  of  the  Markov  time  series  yj  j  are  generated  so  that 


var(Ty,i)  = 


(5) 


This  results  in  a  stationary  time  series. 


A  plot  of  Logio(cty)  against  Logio(x)  for  an  individual  Markov  process  possesses  a  single  maximum  at 
T  =  1.817? .  Plots  of  Logio(cty)  against  Logio(x)  obtained  from  the  sum  of  several  series  of  Markov 
processes  summed  as  described  in  equation  (4)  results  in  a  plot  of  Logio(ay)  against  Logio(x)  that 
exhibits  a  “ripple”  effect;  see  Figure  4.  The  FFM  may  be  constructed  so  that  the  underlying  gradient 
approximates  zero.  It  is  relatively  simple  to  construct  a  linear  combination  of  Markov  processes  such 
that  a  plot  of  Logio(cty)  against  Logio(x)  has  a  gradient  between  -1/2  and  1/2 ,  and  so  models  any 
power-law  noise  process  intermediate  in  characteristics  between  WFM  and  RWFM. 
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Figure  2.  An  example  of  a  simulated  FFM  noise  process  and  its  Markov  components 
(Ml  to  M4),  offset  by  10,  20  30  and  40  units,  respectively,  for  clarity. 
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Figure  3.  Plots  of  Logio(Oy)  against  Logio(x)  obtained  from  both  simulation  and 
theory  for  a  FFM  noise  process  and  its  Markov  components  (Ml  to  M4). 


The  effect  of  increasing  the  “spacing  ratio”  s  =  between  the  Markov  processes  from 

RJ^^  =  Rj  /  8  to  Rj^^  =  R  -  /16  is  shown  in  Figure  4.  The  magnitude  and  period  of  the  “ripple” 

observed  on  the  plots  of  Logio(Oy)  against  Logio(t)  is  significantly  increased.  There  is,  however,  a 
price  to  be  paid  in  using  a  small  spacing  ratio  in  terms  of  additional  computational  time.  The  “ripple” 


284 


55**  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


is  not  clearly  observed  on  the  FFM  curves  of  Figure  3  due  to  the  y-axis  scale  interval  being  much 
larger. 


Figure  4.  Plots  of  Logio(Oy)  against  Logio(x)  obtained  from  theory  when  modelling 
FFM  using  Markov  noise  processes  with  spacing  ratios  of  8  and  16  respectively.  A 
bias  has  been  added  to  the  (ratio  =  16)  curve  for  clarity. 
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Figure  5.  Plots  of  the  autocorrelation  function  p(l)  obtained  as  a  function  of  lag  I 
derived  from  theory  and  by  simulation  for  FFM  noise  processes  and  for  the  four 
Markov  processes  (Ml  to  M4)  used  in  its  construction.  20,000  points  were  used  in 
the  simulation. 
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Figure  5  shows  plots  of  the  autocorrelation  function,  obtained  from  both  theory  and  simulation.  The 
autocorrelation  of  the  FFM  decays  much  slower  than  would  be  expected  from  a  single  exponential 
decay;  this  “long-memory”  characteristic  of  the  FFM  is  well  known  [7].  From  the  above  plots  it  is 
reasonable  to  assume  that  a  good  approximation  to  FFM  has  been  achieved.  It  must  also  be 
remembered  that  it  is  only  possible  to  generate  an  autocorrelation  function  because  we  are  modelling 
FFM  with  a  finite  number  of  Markov  processes;  otherwise,  the  FFM  noise  process  would  not  be 
stationary. 


3.  DESIGN  OF  THE  KALMAN  FILTER 

The  development  of  Kalman  filter  clock  algorithms  have  been  described  previously  [2,4].  In  this 
work,  each  Markov  process  contributes  an  extra  component  to  the  state  vector.  When  using  three 
clocks  where  the  FFM  for  each  clock  is  modelled  by  four  Markov  processes,  the  state  vector  is  given 
by: 


^2,  •••m24,X3,j3,Z3,m3i  •••m34)^ 


(6) 


where  x; ,  yi  and  Zj  are  the  time  offset,  frequency  offset,  and  linear  frequency  drift  offset  between  the 
clock  i  and  an  ideal  clock,  m^  is  the  frequency  offset  between  the  j*  Markov  component  of  the  FFM 
of  clock  i.  The  composite  clock  is  obtained  as  the  Kalman  filter’s  estimate  of  this  “ideal”  clock. 

The  component  (hi;  of  the  state  propagation  matrix  for  clock  i  is  given  by: 


f. 


1 


O  = 

ll 


rl  (l-exp(-RiTo)) 
2 
To 
1 

exp  (-7?,  To) 


(l-exp(-f?4ro))^ 

R, 


(7) 


exp(-7?4ro)  ) 


where  C>  =  Jzag(C>*,..,C)")  ,  there  being  n  clocks  in  the  ensemble.  The  top  left  3-by-3  sub-matrix  is 
well-known  and  the  additional  terms  represent  the  Markov  components  [9]. 


The  component  Qo  of  the  process  covariance  matrix  Q  for  clock  i  is  given  by: 


Qun 

2  7-^/9 
^RWFM^O  '  ^ 


TTrwfmTq  0  Tf^a^2{R^TQ)0'J^^^ 

,2 


0 

Qii  ~  I  ^o^i2(7^iT3)<Tua 


Tla^2{R,T,)ai 


TTrwfmTq  ® 

0  0 


^0^22  (7^1  ^o)*T 


Ma 


70022  (7^4  ^o)*TmC4  J 


(8a) 
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where  Q  =  diag{Q^y,..,  ) ,  (rhffM  =  0  as  we  assume  that  RRFM  is  not  present. 


Qan  ^WFM  ^0  /3  +  ro'a„(i?,ro)cr^ci  +  ..  +  rla,,{R,T,)a^ 


MCA 


«iiW  = 


-  3/2  +  X  +  2  exp  (-x)  -  exp  {-2x)l2 


X 


an  (x)  = 


1/2  -  exp(-x)  +  exp(-2x)/2 


«22W  = 


1  -  exp  (-2x) 
2x 


where  are  the  differential  variances  of  the  Markov  noise  process. 


(8b) 

(8c) 

(8d) 

(8e) 


The  and  terms  of  the  top  left  3-by-3  matrix  are  well-known  [10]  and  the  additional 

terms  describe  the  noise  added  by  the  Markov  processes  [9].  There  is  assumed  to  be  no  correlation 
between  the  noise  observed  in  individual  clocks;,  hence,  the  sub-matrices  Qij  have  only  zero  elements. 


In  this  work,  we  assume  clock  1  to  be  the  reference  clock.  Measurement  are  made  of  (clock  i  -  clock 
1).  The  resulting  design  matrix  in  the  case  of  three  clocks  is  given  by  H,  where: 


H  = 


"-1  0000001000 
^-1  0000000000 


0  0 
0  0 


0  0  0 
0  1  0 


0  0 
0  0 


0  0 
0  0 


0^ 

0. 


(9) 


4.  COVARIANCE  REDUCTION 

The  performance  of  Kalman-filter-based  clock  algorithms  has  been  limited  due  to  the  physical 
parameters  estimated  by  the  state  vector  being  only  partly  observable  from  the  measurements  made 
between  individual  clock  pairs.  Greenhall  [4]  has  recently  overcome  these  problems  through  the 
development  of  a  method  described  as  covariance  reduction.  Following  the  application  of  the  error 
covariance  update  equation,  we  apply  the  following  operation  to  the  parameter  covariance  matrix  P^: 

P^=SP^S^  (10) 


where  the  ith  clock  component  of  S  is  given  by: 

"0 

1 


(11) 
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where  S  =  diagl^S^^  .  S  is  formed  from  an  identity  matrix  where  the  elements  corresponding 

to  the  X  (phase)  component  of  the  state  vector  are  set  to  zero.  The  covariance  reduction  technique  is 
applied  to  all  work  undertaken  in  this  paper. 


5.  EVALUATING  THE  PERFORMANCE  OF  THE  KALMAN 
FILTER  USING  SIMULATED  DATA 

The  performance  of  the  Kalman  filter  for  simulated  data  is  shown  in  Figures  6  to  9  below.  Three 
clocks  are  used  in  each  simulation.  In  Figures  6  and  7,  clocks  1  and  3  exhibit  only  WFM  noise  and 
clock  2  FFM.  In  Figures  8  and  9,  clocks  1  and  3  exhibit  RWFM  and  clock  3  FFM,  the  FFM  being 
modelled  using  four  Markov  processes.  The  clock  phase  outputs  are  shown  in  Figures  6  and  8. 
Figures  7  and  9  show  plots  of  Logio(Cfy)  against  Logio(T). 


Phase  offset 


Figure  6.  Phase  offsets  of  the  three  clocks  used  in  the  simulation.  Clocks  1  and  3 
exhibit  WFM  and  clock  2  FFM.  The  resulting  composite  is  shown. 


Figures  6  to  9  show  the  performance  of  the  individual  clocks  and  the  resulting  composite.  At  almost 
all  averaging  times,  the  stability  of  the  composite  clock  is  equal  to  or  better  than  the  stability  of  the 
most  stable  clock  (Figures  7  and  9).  Data  from  the  first  5%  of  the  time  series  are  not  used  in  the 
computation  of  the  plots  of  Logio(cty)  against  Logio(x)  so  as  to  give  the  Kalman  Filter  time  to  “settle 
down.”  At  any  averaging  time  (x),  it  is  possible  to  form  an  optimally  weighted  mean  of  the  clocks 
that  will  minimize  the  AVAR  of  the  resulting  composite;  these  weights  are  given  by  Wi,  where: 

w  . 

w,  and  W,.  =l/a-^,(r)  (12) 

k=\ 
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Phase  offset  AVAR  plots 


Figure  7.  Plots  of  Logio(Oy)  against  Logio(T)  obtained  from  the  individual  clocks  and 
from  the  composite  obtained  from  the  Kalman  filter.  The  final  two  curves  show  the 
Oy  obtainable  from  an  optimal  weighted  mean  of  the  clock  measurements  and  those 
obtained  from  a  fixed  weighting. 


The  resulting  plot  of  Logio(Cfy)  against  Logio(x)  is  shown  as  the  “optimal”  curve  in  Figures  7  and  9.  It 
is  important  to  note  that,  for  this  curve,  different  clock  weights  have  been  used  at  each  t.  By 
comparing  the  composite  curve  against  the  optimal  curve,  we  conclude  that  the  composite  curve 
approaches  the  optimal  curve  at  all  averaging  times,  demonstrating  that  the  addition  of  Markov 
components  to  the  state  vector  may  result  in  a  clock  algorithm  that  is  close  to  optimal  at  all  averaging 
times.  Finally,  we  take  a  single  clock  weighting  optimal  at  x  =  10,000s.  By  comparing  this  “sample” 
curve  with  the  composite  curve,  we  are  able  to  observe  the  difference  in  performance  between  the 
Kalman-filter-based  algorithm  and  a  simple  weighted  mean  algorithm.  The  Kalman-filter-based  clock 
algorithm  performs  close  to  optimal  at  all  averaging  times,  while  the  simple  weighted  mean  clock 
algorithm  is  optimal  at  a  single  averaging  time,  and  performs  sub-optimally  at  all  other  averaging 
times.  It  must  be  remembered,  however,  that  the  superior  performance  of  the  Kalman-filter-based 
clock  algorithm  may  only  be  obtained  if  the  noise  processes  do  indeed  approximate  a  linear 
combination  of  WFM,  FFM,  and  RWFM  noise  processes  and  requires  that  good  estimates  for  the 
corresponding  noise  parameters  are  available. 


6.  OBSERVING  THE  MARKOV  COMPONENTS  OF  THE  STATE 
VECTOR 

The  successful  operation  of  this  clock  algorithm  is  due  the  Kalman  filter’s  ability  to  correctly  estimate 
the  Markov  components  of  the  state  vector.  Figure  10  shows  the  simulated  data  for  two  of  the 
Markov  processes  M3  and  M4  and  the  resulting  state  vector  estimate.  There  is  a  high  degree  of 
correlation  between  the  two  M4  curves  and,  to  a  lesser  extent,  between  the  two  M3  curves.  These 
results  suggest  that  the  Kalman  filter  is  correctly  estimating  the  Markov  components  of  the  state 
vector. 
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Phase  offset 


Figure  8.  Phase  offsets  of  the  three  clocks  used  in  the  simulation.  Clocks  1  and  3 
exhibit  RWFM  and  clock  2  FFM.  The  resulting  composite  is  shown. 
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Figure  9.  Plots  of  Logio(Oy)  against  Logio(T)  obtained  from  the  individual  clocks  and 
from  the  composite  obtained  from  the  Kalman  filter.  The  final  two  curves  show  the 
Oy  obtainable  from  an  optimal  weighted  mean  of  the  clock  measurements  and  those 
obtained  from  a  fixed  weighting. 
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Markov  noise  simulation  and  state  vector  estimates 
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Figure  10  shows  the  Markov  state  vector  estimates  M3  and  M4  and  corresponding 
simulated  Markov  noise  processes.  Integer  offsets  of  between  zero  and  three  units 
have  been  added  for  clarity. 


7.  ESTIMATING  THE  CLOCK  NOISE  PARAMETERS 


To  operate  this  Kalman  filter  clock  algorithm  using  real  clock  data,  it  is  necessary  to  first  estimate  the 
magnitude  of  the  noise  parameters  ,  and  .  A  method  to  estimate  the  magnitude  of 

noise  parameters  is  described  in  [5].  Estimates  of  a  scaled  AVAR  are  obtained  from  the  measurement 
data  set  at  averaging  times  T  =2" where  «  is  a  positive  integer  and  to  is  the  minimum  spacing  of 

the  data  set.  An  iterative  least  squares  process  is  then  applied  to  produce  a  near  optimal  estimate  of 
the  noise  parameters  in  the  presence  of  highly  correlated  AVAR  estimates. 


In  this  paper,  the  use  of  a  scaled  AVAR  has  been  replaced  with  a  D  (4)  Daubechies  wavelet  variance 
[111.  The  D  (4)  wavelet  is  itself  a  second-difference  variance.  In  this  work,  the  wavelet  variances  are 
applied  to  frequency  measurements,  and  are  used  to  resolve  ,  and  noise 


parameters.  (TppM  represents  a  FFM  noise  process  that  is  approximated  as  the  sum  of  n  Markov 


noise  processes,  each  with  variance 


(T,. 


as  described  in  equation  (2).  The  advantage  of  using 


the  wavelet  variance  is  that,  in  the  presence  of  power-law  noise  processes,  the  correlation  between  the 
estimates  of  individual  variances  is  minimal  and  may  be  neglected.  This  significantly  simplifies  the 
analysis  process,  in  particular  where  there  are  a  relatively  large  number  of  variance  estimates  and 
measurements.  It  must  be  remembered  that  we  are  making  measurements  between  pairs  of  clocks;  the 
noise  parameter  resolution  technique  must  be  combined  with  a  three-cornered  hat.  In  the  case  of  only 
three  clocks  being  present,  we  construct  the  following  least-squares  problem: 
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D  =  H(J  (13b) 

where  Jyk  is  the  D  (4)  wavelet  variance  estimated  from  measurements  of  the  (j-k)  clock  pair.  The 
elements  of  the  design  matrix  H  may  be  simply  computed  for  wavelets  of  small  width  and  obtained 
for  wavelets  of  larger  width  by  extrapolation.  The  variance  of  the  wavelet  variance  var  (dj)  may  be 
estimated  as: 


Yar{d j)  =  2d^.  / n  (14) 

where  n  is  the  number  of  samples  in  the  wavelet  estimate;  these  are  assumed  to  be  uncorrelated. 

An  example  of  the  estimate  of  noise  parameters  is  shown  in  Figure  1 1  below.  The  first  three  curves 
show  plots  of  Logio(ay)  against  Logio(x)  obtained  for  NPL’s  three  active  hydrogen  masers  using  a 
three-cornered  hat.  The  same  data  set  is  then  used  to  estimate  the  ,  <7^^^ ,  and  crlfypM  noise 
parameters  for  each  clock.  These  noise  parameters  are  themselves  used  to  determine  theoretical 
curves  of  Logio(ay)  against  Logio(x);  these  are  shown  in  the  second  three  curves.  The  two  curves 
obtained  for  each  clock  are  in  broad  agreement. 


8.  FIRST  RESULTS  USING  REAL  CLOCK  MEASUREMENTS 

The  Kalman  filter  clock  algorithm  has  been  applied  to  10  days  of  data  from  NPL’s  three  active 
hydrogen  masers.  The  resulting  plots  of  Logio(ay)  against  Logio(x)  obtained  from  (Composite  - 
Individual  Clock)  estimates  are  shown  in  Figure  12.  Because  we  are  unable  to  compare  against  an 
“ideal”  clock,  the  stability  of  (Composite  -  Individual  Clock)  estimates  are  examined.  The  results  are 
encouraging;  however,  it  must  be  remembered  that  there  will  be  a  high  degree  of  correlation  between 
the  most  stable  clocks  used  in  the  algorithm  and  the  composite.  Future  work  will  include  comparing 
the  stability  of  (Composite  -  Individual  Clock)  estimates,  obtained  from  theory  and  those  obtained 
from  applying  the  clock  algorithm  to  real  data. 

In  the  longer  term,  it  is  hoped  to  compare  the  output  of  the  clock  algorithm  against  NPL’s  cesium 
fountain,  as  well  as  against  other  UTC  (k)  time  scales  using  two-way  time  transfer  measurements,  and 
also  to  develop  methods  to  estimate  the  stability  of  the  resulting  composite. 
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PLOTS  OF  LOG^Q(ay)  AGAINST  LOG^(,(t(s)) 


Figure  11.  Plots  of  Logio(Gy)  against  Logio(x)  obtained  for  NPL’s  three  active 
hydrogen  masers  using  a  three-cornered  hat,  and  obtained  using  the  estimated  noise 
parameters. 


Figure  12.  Plots  of  Logio(Oy)  against  Logio(x)  obtained  from  (Composite  -  Individual 
Clock)  estimates  obtained  from  applying  the  Kalman  filter  clock  algorithm  to  NPL’s 
three  active  hydrogen  masers. 
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9  CONCLUSIONS 

We  have  shown  that  Flicker  Frequency  Modulation  noise  may  be  described  approximately  as  a  linear 
combination  of  Markov  noise  processes.  This  description  of  FFM  has  been  successfully  included  in  a 
Kalman-filter-based  clock  algorithm.  The  result  is  a  clock  algorithm  that  will  produce  close  to 
optimal  performance  at  all  averaging  times,  when  operating  in  the  presence  of  FFM.  The  algorithm 
has  been  tested  using  simulated  data.  First  results  have  been  obtained  using  real  clock  data. 
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Questions  and  Answers 

DON  PERCIVAL  (University  of  Washington):  In  statistical  terms,  what  you  are  doing  is  looking  at  a 
linear  combination  of  autoregressive  processes  of  order  one.  Another  approach  which  has  been  taken  in 
the  literature  to  try  to  approximate  flicker  noise  processes  is  to  look  at  higher-order  autoregressive 
processes.  Instead  of  using  nth  first-order  processes,  you  could  think  about  an  nth  fourth-order  auto¬ 
regressive  process.  That  can  be  embedded  in  the  Kalman  filter  very  simply.  Is  there  some  advantage  to 
using  the  linear  combination  of  first-order  processes,  as  opposed  to  the  nth  fourth-order  autoregressive 
process? 

JOHN  DAVIS:  To  be  quite  honest,  it  is  just  the  simplest  way  I  could  see  how  to  do  it,  adding  up  the 
simple  Markov  processes.  We  have  not  looked  at  anything  else.  If  we  took  a  view  after,  I  could  certainly 
try  to  do  that.  That  might  be  interesting. 
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